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Nonlinear d 10 -ML 2 Transition-Metal Complexes 

Lando P. Wolters [a] and F. Matthias Bickelhaupt* [ab] 



We have investigated the molecular geometries of a series of 
dicoordinated d 10 -transition-metal complexes ML 2 (M = Co~, 
RhT, \r, Ni, Pd, Pt, Cu + , Ag + , Au + ; L = NH 3 , PH 3/ CO) using rela- 
tivists density functional theory (DFT) at ZORA-BLYP/TZ2P. Not 
all complexes have the expected linear ligand-metal-ligand 
(L-M-L) angle: this angle varies from 180° to 128.6° as a func- 
tion of the metal as well as the ligands. Our main objective is 
to present a detailed explanation why ML 2 complexes can 
become bent. To this end, we have analyzed the bonding 
mechanism in ML 2 as a function of the L— M— L angle using 



quantitative Kohn-Sham molecular orbital (MO) theory in com- 
bination with an energy decomposition analysis (EDA) scheme. 
The origin of bent L— M— L structures is jt backdonation. In sit- 
uations of strong jt backdonation, smaller angles increase the 
overlap of the ligand's acceptor orbital with a higher-energy 
donor orbital on the metal-ligand fragment, and therefore 
favor jt backdonation, resulting in additional stabilization. The 
angle of the complexes thus depends on the balance between 
this additional stabilization and increased steric repulsion that 
occurs as the complexes are bent. 



Introduction 

Dicoordinated d 10 -transition-metal complexes ML 2 occur in nu- 
merous catalytic reaction mechanisms. [1] These complexes, in 
general, have a linear geometry [2 ' 3,4,5] with a ligand-metal- 
ligand (L-M-L') angle (or bite angle) of 180°, although excep- 
tions^' 71 have been observed. This geometrical preference can 
be easily understood for a closed-shell d 10 configuration. In 
most cases, the dominant bonding orbital interaction is a dona- 
tion from the ligand's lone-pair orbitals into the empty metal 
(n + 1)s atomic orbital (AO), which has a ligand-metal bond 
overlap that is independent of the L— M-L' angle (see 
Figure 1). [8] At the same time, the steric repulsion associated 
with a L-L' overlap between the lone pairs (and other closed 
shells) of the two ligands yields a force that maximizes their 
mutual distance and thus yields the well-known linear L-M-L' 
arrangement. 

The same conclusion is obtained if one uses valence shell 
electron pair repulsion (VSEPR) theory adapted for treating 
transition-metal complexes/ 9, 101 or more sophisticated methods 
based on molecular orbital (MO) theory. Proceeding from the 
latter, one can deduce the preference for linear over bent ML 2 
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Figure 1. a Donation has no preference (left, middle) whereas sterics favor 
linear L-M-L (right). 



complexes from the number of electrons in the valence orbi- 
tals and the dependence of the orbital energies on the geo- 
metrical parameter of interest (here, the L— M— L angle) in 
Walsh diagrams. [8] These diagrams show again that dicoordi- 
nate d 10 -transition-metal complexes, for example, Ag(NH 3 ) 2 + , 
adopt a linear geometry due to the significant destabilization 
of the metal d xz AO by the ligand's lone-pair orbitals in combi- 
nation with steric repulsion between the latter upon bending 
(see below). Nearly all instances with substantial deviations of 
the L— M— L bite angle from linearity are complexes in which 
this distortion is imposed by the structural constraints in bi- 
dentate ligands in which a bridge or scaffold forces the two 
coordinating centers L towards each other. [1b_d] 

In this work, we show that d 10 -ML 2 complexes are not neces- 
sarily linear and may even have a pronounced intrinsic prefer- 
ence to adopt a nonlinear equilibrium geometry. To this end, 
we have investigated the molecular geometries and electronic 
structure of a series of d 10 -ML 2 complexes (M = Co~, Rh", Ir", 
Ni, Pd, Pt, Cu + , Ag + , Au + ; L = NH 3 , PH 3 , CO) using relativistic 
density functional theory (DFT). Simple d 10 -ML 2 complexes are 
found with substantial deviations from linearity, featuring bite 
angles as small as 131° or even less. All that is necessary for 
bent d 10 -ML 2 complexes to occur is sufficiently strong jt back- 
donation. This emerges from our detailed metal-ligand bond- 
ing analyses in the conceptual framework of quantitative MO 
theory contained in Kohn-Sham DFT. The analyses explain the 
phenomenon and provide a tool for rationally tuning the bite 
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angle. Based on our analyses, we can augment the text-book 
Walsh diagram for bending ML 2 complexes involving only 
a donation with an extended Walsh diagram that also includes 
jt backbonding. 

Theoretical Methods 

Computational details: All calculations were carried out using 
the Amsterdam Density Functional (ADF) program developed 
by Baerends and co-workers [11, 12,131 The numerical integration 
was performed using the procedure developed by te Velde 
et al. [14] The molecular orbitals (MOs) were expanded in a large 
uncontracted set of Slater-type orbitals (STOs): TZ2P (no Gaus- 
sian functions are involved). The TZ2P basis set [15] is of triple-^ 
quality for all atoms and has been augmented with two sets of 
polarization functions, that is, 2p and 3d on H, 3d and 4f on C, 
N, 0 and P, 4p and 4f on Co, Ni, Cu, 5p and 4f on Rh, Pd and 
Ag and 6p and 5f on Ir, Pt and Au. An auxiliary set of s, p, d, f 
and g STOs was used to fit the molecular density and to repre- 
sent the Coulomb and exchange potentials accurately in each 
self-consistent field (SCF) cycle. All electrons are included in 
the variational treatment (no frozen-core approximation used). 

Equilibrium structures were obtained by optimizations using 
analytical gradient techniques. [16] Geometries and energies 
were calculated at the BLYP level of the generalized gradient 
approximation (GGA): exchange is described by Slater's Xa po- 
tential/ 171 with nonlocal corrections due to Becke [18] added self- 
consistently, and correlation is treated using the gradient-cor- 
rected functional of Lee, Yang and Parr. [19] Scalar relativistic ef- 
fects were accounted for using the zeroth-order regular ap- 
proximation (ZORA). [20] This approach has been extensively 
tested and was shown to agree well with high-level coupled- 
cluster reference data. [21] Energy minima have been verified 
through vibrational analysis/ 221 All minima were found to have 
zero imaginary frequencies. The PyFrag program was used to 
facilitate the analyses of the bonding mechanism as a function 
of the L-M-L angle. [23] 

Bond energy analysis: The bond energy AE is decomposed 
into the strain energy A£ strain , that is associated with the geo- 
metrical deformation of the frag- 
ments as the bond formation 
takes place, plus the actual inter- 
action energy AE int between the 
deformed fragments [Equa- 
tion (1)]. 



three physically meaningful terms [Eq. (2)] using a quantitative 
energy decomposition scheme developed by Ziegler and 

Rauk. [25] 



AE int (£) = Al/ elstat (£) + AE Pauli (£) + Af ol (C) 



(2) 



The term Al/ e , stat corresponds to the classical electrostatic in- 
teraction between the unperturbed charge distributions p A W + 
p B (r) of the prepared or deformed fragments A and B (see 
below for definition of the fragments) that adopt their posi- 
tions in the overall molecule AB, and is usually attractive. The 
Pauli repulsion term A£ PauN comprises the destabilizing inter- 
actions between occupied orbitals and is responsible for the 
steric repulsion. This repulsion is caused by the fact that two 
electrons with the same spin cannot occupy the same region 
in space. It arises as the energy change associated with the 
transition from the superposition of the unperturbed electron 
densities p A (r)+p B (r) of the geometrically deformed but isolat- 
ed fragments A and B, to the wavefunction x ¥° = NA[ x ¥ a x ¥q\, 
that properly obeys the Pauli principle through explicit anti- 
symmetrization {A operator) and renormalization (N constant) 
of the product of fragment wavefunctions (see Ref. [24] for an 
exhaustive discussion). The orbital interaction AE oi accounts for 
charge transfer (interaction between occupied orbitals on one 
fragment with unoccupied orbitals on the other fragment, in- 
cluding the HOMO-LUMO interactions) and polarization 
(empty-occupied orbital mixing on one fragment due to the 
presence of another fragment). It can be further divided into 
contributions from each irreducible representation r of the 
interacting system [Eq. (3)]. 

M OI (£) = £ r AE£(£) (3) 

Results and Discussion 

Structure and energetics 

Structural and energetic data emerging from our ZORA-BLYP/ 
TZ2P computations are collected in Tables 1-4. Most ML 2 com- 
plexes have a linear L— M-L angle, which leads to either D 3h - 



AE = AE strain + AE in 



(1) 



The interaction energy A£ int (£) 
between two molecular frag- 
ments is analyzed as a function 
of the bite angle £ in the con- 
ceptual framework provided by 
the Kohn-Sham MO method. [24] 
To this end, it is decomposed in 



Table 1. L-M-L angle [°] and linearization energy AE Nn [kcalmol : ] in dicoordinate d 10 -ML 2 complexes. [i 





Group 9 
L-M-L 


A£ lin [b] 




Group 10 
L-M-L 


AF [b] 




Group 11 
L-M-L 


AE lin 


Co(NH 3 ) 2 


180.0 


0 


Ni(NH 3 ) 2 


180.0 


0 


Cu(NH 3 ) 2 + 


180.0 


0 


Co(PH 3 ) 2 


131.8 


6.4 


Ni(PH 3 ) 2 


180.0 


0 


Cu(PH 3 ) 2 + 


180.0 


0 


Co(CO) 2 " 


128.6 


19.9 


Ni(CO) 2 


144.5 


2.1 


Cu(CO) 2 + 


180.0 


0 


Rh(NH 3 ) 2 - 


180.0 


0 


Pd(NH 3 ) 2 


180.0 


0 


Ag(NH 3 ) 2 + 


180.0 


0 


Rh(PH 3 ) 2 - 


141.2 


2.0 


Pd(PH 3 ) 2 


180.0 


0 


Ag(PH 3 ) 2 + 


180.0 


0 


Rh(CO) 2 


130.8 


10.2 


Pd(CO) 2 


155.6 


0.5 


Ag(CO) 2 + 


180.0 


0 


lr(NH 3 ) 2 


180.0 


0 


Pt(NH 3 ) 2 


180.0 


0 


Au(NH 3 ) 2 + 


180.0 


0 


lr(PH 3 ) 2 


144.1 


2.4 


Pt(PH 3 ) 2 


180.0 


0 


Au(PH 3 ) 2 + 


180.0 


0 


lr(CO) 2 


134.2 


13.4 


Pt(CO) 2 


159.0 


0.6 


Au(CO) 2 + 


180.0 


0 



[a] Computed at ZORA-BLYP/TZ2P. 
ometry. 



[b] Relative energy of the linear ML 2 complex relative to its equilibrium ge- 
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Table 2. M-L bond length [A] and BDE [kcalmol '] in monocoordinate d 10 -ML and dicoordinate d 10 -ML 2 com 
plexes. [a] 





M-L 


BDE 




M-L 


BDE 




M-L 


BDE 


C0NH3 [bc] 


1.845 


217.1 


NiNH 3 [c] 


1.827 


77.0 


CuNH 3 + 


1.911 


70.0 


CoPH 3 " [b ' c] 


1.971 


240.6 


NiPH 3 [c] 


1.979 


88.0 


CuPH 3 + 


2.163 


68.7 


CoCO [bc] 


1.630 


280.6 


NiCO [c] 


1.663 


109.3 


CuCO + 


1.833 


50.2 


Co(NH 3 ) 2 [b ' c] 


1.908 


24.0 


Ni(NH 3 ) 2 [c] 


1.888 


36.2 


Cu(NH 3 ) 2 + 


1.919 


61.1 


Co(PH 3 ) 2 [c] 


2.051 


48.2 


Ni(PH 3 ) 2 [c] 


2.108 


36.3 


Cu(PH 3 ) 2 + 


2.232 


48.0 


Co(CO) 2 [c] 


1.715 


76.3 


Ni(CO) 2 [c] 


1.765 


48.6 


Cu(CO) 2 + 


1.882 


45.0 


RhNH 3 [c] 


2.001 


55.5 


PdNH 3 


2.115 


21.6 


AgNH 3 + 


2.212 


48.7 


RhPH 3 [c] 


2.068 


89.9 


PdPH 3 


2.172 


39.4 


AgPH 3 + 


2.415 


47.9 


RhCO [c] 


1.750 


122.0 


PdCO 


1.861 


47.4 


AgCO + 


2.137 


28.4 


Rh(NH 3 ) 2 [c] 


2.089 


22.6 


Pd(NH 3 ) 2 


2.106 


28.6 


Ag(NH 3 ) 2 + 


2.172 


45.2 


Rh(PH 3 ) 2 [c] 


2.196 


38.2 


Pd(PH 3 ) 2 


2.287 


28.6 


Ag(PH 3 ) 2 + 


2.444 


38.1 


Rh(CO) 2 [c] 


1.866 


58.1 


Pd(CO) 2 


1.949 


34.7 


Ag(CO) 2 + 


2.113 


30.7 


lrNH 3 [c] 


1.967 


85.0 


PtNH 3 [c] 


1.981 


50.1 


AuNH 3 + 


2.085 


71.4 


lrPH 3 [c] 


2.056 


126.5 


PtPH 3 [c] 


2.095 


77.3 


AuPH 3 + 


2.240 


84.2 


IrCO [c] 


1.734 


166.3 


PtCO [c] 


1.776 


87.9 


AuCO + 


1.927 


55.0 


lr(NH 3 ) 2 [b ' c] 


2.071 


23.6 


Pt(NH 3 ) 2 [c] 


2.061 


41.6 


Au(NH 3 ) 2 + 


2.088 


64.6 


lr(PH 3 ) 2 [c] 


2.190 


44.1 


Pt(PH 3 ) 2 [c] 


2.249 


38.7 


Au(PH 3 ) 2 + 


2.351 


52.6 


lr(CO) 2 " [c] 


1.854 


66.3 


Pt(CO) 2 [c] 


1.911 


47.1 


Au(CO) 2 + 


2.002 


40.4 



[a] Computed at ZORA-BLYP/TZ2P. Bond dissociation energies (BDEs) are given for the complexes in the elec 
tronic configuration corresponding to a d 10 s° electron configuration and relative to closed-shell d 10 s° meta 
atoms, [b] The d 10 s°-type configuration is an excited state of the complex, [c] The d 10 s° configuration is an excit- 
ed state of the atomic metal fragment. 



symmetric complexes M(NH 3 ) 2 and M(PH 3 ) 2 or D^-sym metric 
complexes M(CO) 2 . However, numerous significantly smaller 
angles appear throughout Table 1 as well, where the symmetry 
of the complexes is lowered to C 2v . For instance, the complexes 
become increasingly bent when the ligands are varied along 
NH 3 (a strong a donor), PH 3 (a a donor and jt acceptor) and 
CO (a strong jt acceptor). This is most clearly seen for the 
group 9 complexes, where, for example, the angle decreases 
along Rh(NH 3 ) 2 ", Rh(PH 3 ) 2 - and Rh(CO) 2 " from 180.0° to 141.2° 
and 130.8° (Figure 2). In a later section, we will show that the 
jt-backbonding properties of the complexes constitute a promi- 
nent part of the explanation of why d 10 -ML 2 complexes can 
adopt nonlinear geometries. The increasingly strong jt back- 
bonding along this series also results in stronger metal-ligand 
bonds (see Table 2 for bond dissociation energies (BDEs) and 
Table 3 for energy decomposition analyses (EDA) results for ML 
complexes). 

The extent of bending systematically decreases when the jt- 
backbonding capability of the metal center decreases from the 
group 9 anions, via neutral group 10 atoms, to the group 11 
cations. This is clearly displayed by the series of isoelectronic 
complexes Rh(CO) 2 ", Pd(CO) 2 and Ag(CO) 2 + along which the 
L-M-L angle increases from 130.8° to 155.6° to 180° (Table 1). 




180° 




141 




131 



Figure 2. Equilibrium geometries computed at ZORA-BLYP/TZ2P. From left to right: 
Rh(NH 3 ) 2 -, Rh(PH 3 ) 2 and Rh(CO) 2 . 



The data in Table 3 for the corre- 
sponding monocoordinate 
RhCO", PdCO and AgCO + nicely 
show how along this series the 
distortive Jt-orbital interactions 
AE*. indeed become weaker, 
from -120 to -51 to -11 kcal 
mol" 1 , respectively. In the case 
of group 9 metals, both phos- 
phine and carbonyl complexes 
are bent, whereas, for group 10 
metals, only the carbonyl com- 
plexes deviate from linearity. 
Complexes with a metal center 
from group 11 all have a linear 
L— M— L configuration. The re- 
duced jt backbonding also leads 
to weaker metal-ligand bonds. 
For the cationic metal centers, 
for which jt backdonation plays 
a much smaller role, the metal- 
ligand BDEs decrease in the 
order NH 3 >PH 3 >CO (see 
Table 2). This trend originates di- 
rectly from the o-donating capa- 
bilities of the ligands as reflected 
by the energy of the lone-pair orbital e(LP), which decreases in 
this order (see Table 4). Note that, for the same reason, the ba- 
sicity of the ligand as measured by the proton affinity (PA) de- 
creases along NH 3 > PH 3 > CO. [26] For the anionic group 9 metal 
centers, the opposite order is found, that is, metal-ligand BDEs 
decrease in the order CO> PH 3 > NH 3 , following the jr-accept- 
ing capabilities of the ligands. 

Linearity also increases if one descends in a group. For ex- 
ample, from Ni(CO) 2 to Pd(CO) 2 to Pt(CO) 2 , the L-M-L angle in- 
creases from 144.5° to 155.6° to 159.0°. Interestingly, this last 
trend is opposite to what one would expect proceeding from 
a steric model. If one goes from a larger to a smaller metal 
center, that is, going up in a group, the ligands are closer to 
each other and thus experience stronger mutual steric repul- 
sion. But instead of becoming more linear to avoid such repul- 
sion, the complexes bend even further in the case of the small- 
er metal. For example, when the palladium atom in Pd(CO) 2 is 
replaced by a smaller nickel atom, the L-M-L angle decreases 
from 155.6° in Pd(CO) 2 to 144.5° in Ni(CO) 2 . Later on, we show 
that this seemingly counterintuitive trend also originates from 
enhanced jt backbonding which dominates the increased 
steric repulsion. 

General bonding mechanism 

The bending of our model complexes can be under- 
stood in terms of a monocoordinate complex to 
which a second ligand is added either in a linear or 
a bent arrangement, ML + L^ML 2 (see below). Using 
Pd(CO) 2 as an example, we start from a PdCO frag- 
ment, and consider the addition of the second CO 
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Table 3. Energy decomposition analyses [kcalmol 1 ] and orbital energies e [eV] for the metal-ligand bonds in 
monoligated transition-metal complexes M-L [a] 



ML 


Af 


A£ int 


Al/elstat 


Afpauli 




AE° 0] 


^£7i [b] 


c[dj 




c[dj 


C0NH3 
C0PH3- 
CoCCT 


-217.1 
-240.6 
-280.6 


-218.4 
-241.7 
-286.4 


-110.0 
-197.9 
-233.4 


166.3 
204.5 
274.5 


-274.7 
-248.2 
-327.5 


-241.8 
-123.9 
-141.7 


-32.9 
-124.4 
-185.8 


+ 1.84 
+ 1.67 
+ 1.34 


+ 2.91 
+ 1.81 
+ 1.17 


+ 3.99 
+ 3.38 
+ 3.20 


RhNH 3 
RhPH 3 
RhCO 


-55.5 
-89.9 
-122.0 


-56.2 
-90.3 
-126.0 


-143.2 
-269.7 
-273.3 


202.1 
311.7 
364.1 


-115.1 
-132.3 
-216.8 


-110.8 
-61.7 
-96.7 


-4.3 
-70.6 
-120.1 


+ 1.72 
+ 1.49 
+ 1.05 


+ 1.83 
+ 0.91 
-0.09 


+ 2.53 
+ 2.20 
+ 1.56 


lrNH 3 
lrPH 3 
IrCO 


-85.0 
-126.5 
-166.3 


-85.8 
-127.2 
-171.3 


-196.9 
-349.2 
-353.5 


268.9 
396.0 
461.5 


-157.8 
-174.1 
-279.2 


-142.9 
-85.9 
-129.6 


-14.9 
-88.2 
-149.7 


+ 1.54 
+ 1.18 
+ 0.63 


+ 2.16 
+ 0.73 
-0.26 


+ 2.91 
+ 2.28 
+ 1.68 


NiNH 3 
NiPH 3 
NiCO 


-77.0 
-88.0 
-109.3 


-77.3 
-88.7 
-110.4 


-116.2 
-161.3 
-171.6 


139.8 
173.3 
210.3 


-100.8 
-100.7 
-149.1 


-94.5 
-50.8 
-60.4 


-6.3 
-49.9 
-88.7 


-3.28 
-3.79 
-4.89 


-2.99 
-3.93 
-5.40 


-2.21 
-2.90 
-4.14 


PdNH 3 
PdPH 3 
PdCO 


-21.6 
-39.4 
-47.4 


-21.7 
-39.8 
-47.8 


-88.0 
-166.2 
-161.4 


105.1 
190.3 
213.3 


-38.8 
-63.8 
-99.7 


-34.5 
-35.3 
-48.0 


-4.4 
-28.5 
-51.8 


-3.46 
-4.49 
-5.28 


-3.81 
-5.29 
-6.48 


-3.47 
-4.56 
-5.53 


PtNH 3 
PtPH 3 
PtCO 


-50.1 
-77.3 
-87.9 


-50.4 
-78.9 
-88.7 


-170.1 
-273.9 
-271.6 


211.4 
310.3 
356.9 


-91.7 
-115.3 
-174.0 


-82.0 
-70.5 
-91.6 


-9.7 
-44.8 
-82.4 


-4.19 
-4.92 
-5.97 


-4.46 
-5.72 
-7.28 


-3.72 
-4.53 
-5.77 


CuNH 3 + 
CuPH 3 + 
CuCO + 


-70.0 
-68.7 
-50.2 


-70.1 
-73.5 
-50.3 


-104.5 
-101.7 
-89.8 


86.0 
94.0 
100.7 


-51.7 
-65.8 
-61.2 


-41.9 
-51.8 
-38.8 


-9.8 
-14.0 
-22.4 


-11.80 
-11.99 
-13.7 


-12.13 
-12.44 
-14.28 


-12.02 
-12.15 
-13.90 


AgNH 3 + 
AgPH 3 + 
AgCO + 


-48.7 
-47.9 
-28.4 


-48.7 
-51.8 
-28.6 


-73.3 
-84.3 
-59.1 


58.8 
81.3 
67.2 


-34.2 
-48.8 
-36.7 


-28.5 
-39.9 
-26.2 


-5.8 
-8.9 
-10.6 


-12.56 
-12.41 
-14.08 


-13.60 
-13.67 
-15.07 


-13.57 
-13.85 
-14.86 


AuNH 3 + 
AuPH 3 + 
AuCO + 


-71.4 
-84.2 
-55.0 


-71.6 
-91.0 
-55.1 


-124.8 
-177.9 
-149.0 


123.2 
187.2 
188.4 


-70.0 
-100.3 
-94.5 


-60.3 
-80.9 
-64.9 


-9.7 
-19.4 
-29.7 


-12.49 
-12.52 
-14.20 


-13.32 
-13.70 
-15.53 


-12.92 
-13.06 
-14.73 



[a] Computed at ZORA-BLYP/TZ2P. See [Eqs. (1)-(3)]. [b] Also includes small contributions from 6 orbital inter- 
actions, which can only be separated for C^-sym metric MCO complexes. There, the 6 term amounts at most 
to 3.5% of the jt term. 



Table 4. Ligand orbital energies 


e [eV] and proton affinities [kcalmol ] ]. [a] 


c(LP) 


s{n*) PA 


NH 3 -6.05 
PH 3 -6.63 
CO -8.93 


+ 1.42 +201.4 
-0.24 +185.2 
-1.92 +141.5 


[a] Computed at ZORA-BLYP/TZ2P. LP: lone pair, jt*: acceptor orbital. 
Proton affinities (PA) from enthalpies at 298.15 K and 1 atm. 



tailed energy decomposition 
consolidate this picture. 



with the lone pair on the ligand, 
resulting in a "d 0 " orbital that is 
relatively high in energy, at 
-5.3 eV. 

When the second CO ligand 
coordinates opposite the first 
one (i.e., in a linear L-M-L ar- 
rangement), its jt*-acceptor orbi- 
tals interact with the orbitals 
on the PdCO fragment. The 
latter are already considerably 
stabilized by jt backdonation to 
the first CO ligand (Figure 3 B, 
left). When, instead, the second 
ligand is added at an angle of 
90°, its jt* orbitals overlap with 
only one d^ orbital, and with 
one d§ orbital (Figure 3 B, right). 
This d 8 orbital is essentially 
a pure metal d orbital that has 
not yet been stabilized by any 
coordination bond. Consequent- 
ly, this orbital has a higher 
energy and is, therefore, a more 
capable donor orbital for jt back- 
donation into the jt* orbital of 
the second CO ligand. This re- 
sults in a stronger, more stabiliz- 
ing donor-acceptor interaction 
of this pair of orbitals in the 90° 
(Figure 3 A, right) than in the 
180° ML 2 geometry (Figure 3 A, 
left: cf. red-highlighted Jt inter- 
actions). o-Donation interactions 
are affected less by bending. It is 
therefore jt backdonation that 
favors bending. The more de- 
analyses in the following sections 



ligand both at a 180° angle and a 90° angle. Our Kohn-Sham 
MO analyses show that, in PdCO, the degeneracy of the five 
occupied d orbitals on palladium is lowered by interactions 
with the ligand (see Figure 3). Choosing the M— L bond along 
the zaxis, the d xz and d yz orbitals act as donor orbitals for jt 
backdonation into the two jr*-acceptor orbitals on the CO 
ligand, resulting in two stabilized "d K " orbitals at —6.5 eV 
(value not shown in Figure 3). The d xy and d x2 _ y2 (or "d 8 ") orbi- 
tals at -5.5 eV do not overlap and interact with the ligand. 
The d Z 2 orbital is destabilized due to the antibonding overlap 



Bonding mechanism: Variation of ligands 

To understand the trends in nonlinearity of our ML 2 complexes 
(see above and Table 1), we have quantitatively analyzed the 
metal-ligand bonding between ML and the second ligand L as 
a function of the L— M-L angle. The results are collected in 
Table 2 and displayed in Figure 4-7. Most of our model com- 
plexes have a d 10 -type ground-state configuration but not all 
of them, as indicated in detail in Table 2. Yet, all model systems 
discussed here have been kept in d 10 -configuration, to achieve 
a consistent comparison and because, on the longer term, we 
are interested in understanding more realistic dicoordinated 
d 10 -transition-metal complexes that feature, for example, as 
catalytically active species in metal-mediated bond activation. 
We start in all cases from the optimal linear ML 2 structure (i.e., 
the complex optimized in either or D 3h symmetry) and 
then analyze the bonding between ML and L' as a function of 
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Figure 3. A) Schematic MO diagrams for the bonding mechanism between PdCO and CO in linear Pd(CO) 2 (left) 

and at a L-M-L angle of 90° (right): dominant interactions ( ), other interactions ( ), jt backbonding ( ). 

B) Schematic representation of the bonding overlaps of the donating orbital on PdCO (black) with the jr-accepting 
orbital on the second CO ligand (red). 



the L-M-L angle, from 180° to 90°, while keeping all other ge- 
ometry parameters frozen. The analyses were done in C s sym- 
metry, bending the complexes in the mirror plane, with the 
out-of-plane hydrogen atoms of M(NH 3 ) 2 and M(PH 3 ) 2 towards 
each other. Thus, we are able to separate the orbital interac- 
tions symmetric to the mirror plane (A' irrep) from the orbital 
interactions asymmetric to the mirror plane (A" irrep): AE oi (£) = 
A£*(£) + A£^(£) [Eq. (3)]. The use of frozen fragment geome- 
tries allows us to study purely how the interaction energy 
changes as the angle is varied, without any perturbation due 
to geometrical relaxation. Therefore, any change in AE stems 
exclusively from a change in AE int = Al/ e , stat + AE PauM + AE*j + 
AEqi". Note that rigid bending of the linearly optimized L— M— L 
complexes causes minima on the energy profiles to shift to 
larger angles than in fully optimized complexes, but this does 
not alter any relative structural or energy order. 



In Figure 4, we show the 
energy decomposition analyses 
[Eq. (2)] and how they vary 
along the palladium complexes 
Pd(NH 3 ) 2 , Pd(PH 3 ) 2 and Pd(CO) 2 . 
Upon bending the LM— L' com- 
plex from 180° to 90°, the aver- 
age distance between the elec- 
tron density on LM and the 
nuclei of L' decreases (the Pd— P 
distance however remains con- 
stant), which results in a more 
stabilizing electrostatic attraction 
AVeistat- Likewise, the Pauli repul- 
sion AE PauM increases because of 
a larger overlap of the lone pair 
on L' with the d z2 -derived d Q or- 
bital on the ML fragment. The 
latter is the antibonding combi- 
nation of the metal d z2 orbital 
and the ligand lone pair, with 
a fair amount of metal s character 
admixed in an L— M bonding 
fashion. The resulting hybrid or- 
bital is essentially the d z2 orbital 
with a relatively large torus. The 
increase in Pauli repulsion that 
occurs as the L-M-L' angle de- 
creases stems largely from the 
overlap of the lone pair on the 
second ligand L' with this torus. 
For Pd(CO) 2 for example, the 
overlap of the L' lone pair with 
the d a hybrid orbital on ML in- 
creases from 0.05 to 0.28 upon 
bending from 180° to 90°. We 
note that this repulsion induces 
a secondary relaxation, showing 
up as a stabilizing AE^, by 
which it is largely canceled 
again. The mechanism through 
which this relief of Pauli repulsion happens is that, in the anti- 
bonding combination with the L' lone pair, the d Q orbital is ef- 
fectively pushed up in energy and (through its L'-lone-pair 
component) interacts in a stabilizing fashion with the metal 
s-derived LUMO on ML. 

The aforementioned jt backbonding that favors bending 
(see Figure 3) shows up in an increased stabilization in the 
asymmetric AE^" component as the L-M— L angle decreases. To 
more clearly reveal the role of the orbital interactions within 
A" symmetry, we separate the interaction energy AE int into the 
corresponding term AEj- plus the remaining interaction 
energy AE' int , which combines the other interaction terms com- 
prising electrostatic attraction AV ehtaV Pauli repulsion AE PauN , 
and the symmetric orbital interactions AE^: 
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Figure 4. Energy decomposition analysis [Eq. (2)] of the interaction between 
PdL and L in dicoordinated palladium complexes Pdl_ 2 as a function of the 
L-M-L angle (L = NH 3 , PH 3 , CO). 



AE int (£)= AV etat (£) + AE Pauli (C) + A£*(£) + AE*" 

Thus, the interaction energy is split into two contributions 
which are both stabilizing along a large part of the energy pro- 
files studied and which vary over a significantly smaller range. 
Therefore, this decomposition allows us to directly compare 
the importance of AE% with respect to the combined influence 
of all other terms, contained in AF int . The latter contains the 
aforementioned counteracting and largely canceling terms of 
strong Pauli repulsion between A' orbitals and the resulting 
stabilizing relaxation effect AE*. 

The results of this alternative decomposition appear in 
Figure 5, again for the series of palladium complexes Pd(NH 3 ) 2 , 
Pd(PH 3 ) 2 and Pd(CO) 2 . In each of these complexes, bending 
begins at a certain point to weaken the AF int energy term and, 
at smaller L-M-L angles, makes it eventually repulsive as the 
Pauli repulsion term becomes dominant (see also Figure 4). 
Numerical experiments, in which we consider the rigid bend- 
ing process of a complex in which the metal is removed, show 
that steric repulsion between ligands does contribute to this 
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Figure 5. Energy decomposition analysis [Eq. (4)] of the interaction between 
PdL and L in dicoordinated palladium complexes PdL 2 as a function of the 
L-M-L angle (L = NH 3 , PH 3 , CO). 



repulsion, especially at smaller angles. Thus, direct Pauli repul- 
sion between L and L' in LM— L' goes, upon bending from 180° 
to 90°, from 0.3 to 4.6 kcalmol" 1 for Pd(NH 3 ) 2 and from 0.4 to 
9.0 kcalmol" 1 for Pd(CO) 2 (data not shown in Figures). This 
finding confirms that ligands avoid each other for steric rea- 
sons, but it also shows that the effect is small as compared to 
the overall change in the AE int curves (see Figure 5). The domi- 
nant term that causes AE int to go up in energy upon bending 
is the increasing Pauli repulsion that occurs as the L' lone pair 
overlaps more effectively with the LM d G orbital. 

In a number of cases, the stabilization upon bending from 
the asymmetric orbital interactions AE^" dominates the desta- 
bilization from the AF int term. These cases are the complexes 
that adopt nonlinear equilibrium geometries. This AE% term 
gains stabilization upon bending LM-L' because the ^-accept- 
or orbital on the ligand L' moves from a position in which it 
can overlap with a ligand-stabilized LM d^ orbital to a more or 
less pure metal and, thus, up to 1 eV higher-energy d 8 orbital 
(see Table 3), which leads to a more stabilizing donor-acceptor 
orbital interaction (see Figure 5). The gain in stabilization of 
AE^" upon bending and, thus, the tendency to bend increases 
along NH 3 to PH 3 to CO. The reason is the increasing jt-accept- 
ing ability of the ligands as reflected by the energy e(jt*) of 
the ligands' jt* orbital which is lowered from +1.42 to —0.24 
to -1.92 eV, respectively (see Table 4). Thus, for Pd(NH 3 ) 2 , 
where jt backdonation plays essentially no role, the AE^' term 
is stabilized by less than 0.5 kcalmol" 1 if we go from 180° to 
90°. For PH 3 , known as a moderate Jt-accepting ligand, this 
energy term is stabilized by 1.5 kcalmol" 1 from 180° to 90° 
and, for CO, this stabilization amounts to 2.5 kcalmol" 1 . Thus, 
in the case of palladium complexes, the energy profile for 
bending the complexes becomes progressively more flat as 
the ligands are better jt acceptors, but only the carbonyl 
ligand generates sufficient stabilization through increased jt- 
backbonding in AE*." to shift the equilibrium geometry to an 
angle smaller than 180°. 

Bonding mechanism: Variation of metals 

Applying the same analysis along the series Rh(CO) 2 ", Pd(CO) 2 
and Ag(CO) 2 + , reveals a similar but clearer picture (Figure 6). 
Along this series of isoelectronic complexes, the equilibrium 
geometries have L-M-L angles of 130.8°, 155.6° and 180.0°. 
Similar to the results obtained for the series discussed above, 
we again find a AF int term that is relatively shallow and even- 
tually, at small angles, dominated by the Pauli repulsion. The 
A£' int term does not provide additional stabilization upon 
bending the complex. We do observe, however, a AE% compo- 
nent that, from Rh(CO) 2 " to Pd(CO) 2 to Ag(CO) 2 + , becomes 
more stabilizing and also gains more stabilization upon bend- 
ing from 180° to 90°. That is, whereas for Ag(CO) 2 + the AE% 
remains constant at a value of —5.4 kcalmol" 1 as the complex 
is bent from 180° to 90°; the same component for Pd(CO) 2 
starts already at a more stabilizing value of —15.1 kcalmol" 1 at 
180° and is stabilized more than 2.5 kcalmol" 1 as the complex 
is bent to 90°. For Rh(CO) 2 ", the effect of the additional stabili- 
zation upon bending is strongest, almost 10 kcalmol" 1 , as AE^" 
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Figure 6. Energy decomposition analysis [Eq. (4)] of the interaction between 
MCO and CO in dicarbonyl-transition-metal complexes M(CO) 2 as a function 
of the L-M-L angle (M = Rh , Pd, Ag + ). 



Figure 7. Energy decomposition analysis [Eq. (4)] of the interaction between 
MCO and CO in dicarbonyl-transition-metal complexes M(CO) 2 as a function 
of the L-M-L angle (M = Ni, Pd, Pt). 



goes from -28.4 kcalmol" 1 at 180° to -37.3 kcalmol" 1 at 90°. 
The mechanism behind this trend is that the donor capability 
of the metal d orbitals increases as they are pushed up in 
energy from the cationic AgCO + to the neutral PdCO to the 
negative RhCO" (see Table 3). This trend of increasing d orbital 
energies leads to a concomitant strengthening jt backdonation 
and, thus, an increasing energy difference in the LM fragment 
between the pure metal d§ and the ligand-stabilized d^ orbi- 
tals. Thus, the "fresh" d 8 orbitals are higher in energy than the 
ligand-stabilized d,, orbitals by 0.21 to 0.96 to 1.65 eV along 
AgCO + , PdCO and RhCO", respectively (see Table 3). Conse- 
quently, the LM-L/ complexes benefit progressively along this 
series from increasing the overlap of L/ jt* with the higher- 
energy d 8 orbitals in the bent geometry. 

Variation of the metal down a group goes with a less pro- 
nounced increase of the L-M-L angle that originates from 
more subtle changes in the bonding mechanism. The largest 
variation in bite angle is observed along the group 10 com- 
plexes Ni(CO) 2 , Pd(CO) 2 and Pt(CO) 2 which show L-M-L angles 
of 144.5°, 155.6° and 159.0°, respectively (see Table 1). Two fac- 
tors lie behind this trend: (1) a weakening in jt backbonding as 
the metal orbital energy decreases from nickel 3d to palladium 
4d; (2) a steeper increase upon bending in Pauli repulsion be- 
tween PtCO d Q (that has a large torus due to strong admixture 
of the relativistically stabilized Pt 6s AO) and the lone pair of 
the other CO ligand. As shown in Figure 7, the Jt-backbonding 
stabilization of AE^ upon bending is indeed stronger for 
Ni(CO) 2 than for Pd(CO) 2 and Pt(CO) 2 . The difference between 
the latter is small because the greater (more favorable) overlap 
of the jt* orbitals on the ligand with the more extended plati- 
num d orbitals on PtCO compensates for the lower (less favor- 
able) platinum d orbital energy. Figure 7 also shows how the 
AF int term containing the aforementioned Pauli repulsion be- 
comes more rapidly destabilizing at smaller angles for Pt(CO) 2 
than for Ni(CO) 2 and Pd(CO) 2 . Likewise, in the case of group 9 
complexes, the more steeply increasing Pauli repulsion of the 
ligand lone pair with the large iridium d 0 torus pushes the 
equilibrium L-M-L angle of lr(CO) 2 " (134.2°) to a larger value 
than for Rh(CO) 2 " (130.8°; see Table 1). Interestingly, here, the 



linearization energy AE Mn is nevertheless higher for the less 
bent lr(CO) 2 - (13.4 kcal mol" 1 ) than for Rh(CO) 2 " (10.2 kcal 
mol" 1 ) because of the more favorable jt-backbonding overlap 
between IrCO" and CO (see Table 1). This illustrates the subtle- 
ty of the interplay between the two features in the bonding 
mechanism. 

Walsh diagrams 

Based on detailed Kohn-Sham MO analyses of individual com- 
plexes, we have constructed generalized Walsh diagrams corre- 
sponding to bending the ML 2 complexes from 180° to 90°. 
This choice comes down to an alternative perspective on the 
same problem, and the emerging electronic mechanism, why 
bending may occur, is fully equivalent to the one obtained in 
the above analyses based on two interacting fragments 
LM + L', namely: Bending ML 2 to a nonlinear geometry enables 
ligand jt* orbitals (if they are available on L) to overlap with 
and stabilize metal d orbitals that are not stabilized in the 
linear arrangement. The spectrum of different bonding situa- 
tions has been summarized in two simplified diagrams that 
correspond to two extreme situations: weakly jt-accepting li- 
gands (Figure 8A) and strongly jt-accepting ligands (Fig- 
ure 8B). In these diagrams, we position the d z2 orbital in linear 
ML 2 above the other d orbitals, a situation that occurs, for ex- 
ample, for Pd(PH 3 ) 2 . The relative position of the d z2 may 
change, and in some complexes, such as, Rh(NH 3 ) 2 ", it is locat- 
ed below the other d orbitals. These variations do not affect 
the essential property of the orbitals, namely, their change in 
energy upon bending the ML 2 complex. Furthermore, we 
speak about weakly jt-accepting ligands, not just about 
(purely) o-donating ligands, because it turns out that none of 
our model ligands has negligible Jt-accepting capability. The 
resulting Walsh diagrams summarize our results in a more easy 
to use pictorial manner which, in particular for the situation 
with strongly jt-accepting ligands, is novel. 

We first examine the diagram with weakly jt-accepting li- 
gands (Figure 8A). Bending ML 2 from linear to nonlinear signif- 
icantly destabilizes the d xz orbital because of turning on over- 
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Figure 8. Simplified Walsh diagrams for bending ML 2 complexes that emerge from our Kohn-Sham MO analyses (+/- indicate bonding/antibonding) A) with- 
out and B) with jt backbonding. 



lap with the out-of-phase combination of ligand lone pairs. 
This effect is related to the overlap between the LM d a torus 
and the L' lone pair in the fragment approach (see above). At 
small angles, direct ligand-ligand antibonding becomes impor- 
tant. The d Z 2 orbital is slightly stabilized in the nonlinear situa- 
tion due to a decreasing antibonding overlap with the in- 
phase combination of ligand lone pairs, augmented by admix- 
ing with the d X 2_ y2 orbital (see a detailed scheme of this inter- 
mixing in Figure S1 of the Supporting Information). Note that 
if our model ligands would have been purely a donating, the 
d xz , d yz and d xy levels would not be affected by L-M-L bend- 
ing. Yet, they are, although only slightly so. This is a manifesta- 
tion of some Jt backbonding, which is discussed in more detail 
below for the strongly jt-accepting ligands. 

In the case of strongly jt-accepting orbitals (Figure 8 B), 
bending ML 2 from linear to nonlinear still goes with significant 
destabilization of d xz and slight stabilization of d z2 (for the 
same reasons as discussed above for weakly jt-accepting li- 
gands). Jt Backbonding stabilizes both d xz and d yz in the linear 
L-M— L arrangement; bending reduces Jt overlap which causes 
also d yz to increase in energy. A striking phenomenon in the 
ML 2 Walsh diagram with strongly jt-accepting ligands is the 
significant stabilization of the d x2 _ y2 and d xy orbitals that occurs 
as bending moves ligand Jt* orbitals in the right orientation 
for jt-accepting overlap with these orbitals. The resulting stabi- 
lization, if strong enough, can overcome the destabilization of 
the d xz orbital and accounts for the observed bent complexes 
described in this work. This effect is related to the overlap be- 
tween the LM d 8 orbital and the L' jt* in the fragment ap- 
proach (see above). The same effect also nicely accounts for 
the nonlinear structures observed in earlier studies for d° metal 
complexes with n-donating ligandsJ 27,28,29,30,311 For these com- 
plexes, a Jt-bonding mechanism has been proposed in which 
bending is favorable because it effectively increases the 
number of d orbitals that have non-zero overlap with the 
jt-donating orbitals on the ligands. [28] 



Conclusion 

Dicoordinated d 10 -transition-metal complexes ML 2 can very 
well adopt nonlinear geometries with bite angles that deviate 
significantly from the usually expected 180°. This follows from 
our relativistic density functional theory (DFT) computations 
on a broad range of archetypal d 10 -ML 2 model systems. The 
smallest bite angle encountered in our exploration among 
27 model systems amounts to 128.6° for Co(CO) 2 ". 

Nonlinear geometries appear to be a direct consequence of 
jt backbonding. The geometry of d 10 -ML 2 complexes results 
from two opposing features in the bonding mechanism, which 
we have analyzed in terms of the interaction between ML and 
L as a function of the L— M— L angle using quantitative molecu- 
lar orbital (MO) theory and energy decomposition analyses: 
Bending destabilizes the interaction AE int between ML and L 
through increasing steric (Pauli) repulsion between the ligands' 
lone-pair orbital lobes as well as a destabilization, by the latter, 
of the ML d a hybrid orbital; however, bending can also stabilize 
A£j nt because of enhanced Jt backdonation. The reason is that 
the Jt-accepor orbital on the ligand L (e.g., CO jt*) interacts in 
the linear arrangement with an already stabilized ML d n hybrid 
orbital, whereas in the bent geometry, it enters into a more fa- 
vorable donor-acceptor orbital interaction with an unstabi- 
lized, that is, higher-energy metal d 8 orbital. 

Our analyses complement the existing text-book Walsh dia- 
gram for bending ML 2 complexes [8] with a variant that includes 
metal-ligand Jt backbonding. Our findings also contribute to 
a more rational design of catalytically active and selective ML 2 
complexes. [1,32] 
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